Mon. Not. R. Astron. Soc. 000,[T]|5l(2004) Printed 2 February 2008 (MN IATpX style file v2.2) 



Gamma-ray bursts from synchrotron self-Compton emission 



(N 
> 



Boris E. Stern 1 ' 2,3 * and Juri Poutanen 1 * j 

1 Astronomy Division, P.O.Box 3000, 90014 University ofOulu, Finland 

2 Institute for Nuclear Research, Russian Academy of Sciences, 7a, Prospect 60-letija Oktjabrja, Moscow 117312, Russia 
3 Astro Space Center ofLebedev Physical Institute, Profsoyuznaya 84/32, Moscow 117997, Russia 



Accepted, Received 



o 
o 

<N 
G : 

ABSTRACT 

The emission mechanism of the gamma-ray bursts (GRBs) is still a matter of debates. The 
standard synchrotron energy spectrum of cooling electrons Fg oc _E -1 / 2 is much too soft to 
account for the majority of the observed spectral slopes. An alternative in the form of quasi- 
thermal Comptonization in a high compactness source has difficulties in reproducing the peak 
of the observed photon distribution below a few hundred keV. We show here that for typical 
qq ■ parameters expected in the GRB ejecta the observed spectra in the 20-1000 keV BATSE 

qq | energy range can be produced by inverse Compton scattering of the synchrotron radiation in 

■^J- , a partially self-absorbed regime. If the particles are continuously accelerated/heated over the 

' life-time of a source rather than being instantly injected, a prominent peak develops in their 

distribution at a Lorentz factor 7 ~ 30 — 100, where synchrotron and inverse-Compton losses 
are balanced by acceleration and heating due to synchrotron self-absorption. The synchrotron 
peak should be observed at 10-100 eV, while the self-absorbed low-energy tail with Fe oc E 2 
^ (— I can produce the prompt optical emission (like in the case of GRB 990123). The first Compton 

O ^ scattering radiation by nearly monoenergetic electrons can then be as hard as Fe oc E 1 

reproducing the hardness of most of the observed GRB spectra. The second Compton peak 
should be observed in the high energy gamma-ray band, possibly being responsible for the 
emission detected by EGRET in GRB 941017. A significant electron-positron pair production 
reduces the available energy per particle, moving the spectral peaks to lower energies as the 
burst progresses. The regime is very robust, operates in a broad range of parameter space and 
m ' can explain most of the observed GRB spectra and their temporal evolution. 

Key words: gamma-rays: bursts - gamma-rays: theory - methods: numerical - radiation 
mechanisms: non-thermal - scattering 
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1 INTRODUCTION 

Spectra of the prompt soft gamma-ray emission of gamma-ray 
bursts (GRBs) are still not explained and seem mysterious de- 
spite large theoretical efforts devoted to this problem. Already in 
the 1980-ies it was recognized that synchrotron emission from 
the electrons injected at high energies produces cooling spec- 
tra Fe oc E' 1 / 2 (described by a photon spectral index a = 
—3/2) which are muc h too soft to be consistent with that ob - 
served from GRBs (e.g. Bussard 1984; Imamura & Epstein 1987). 
The problem became acute when Pre ece et alJ j2000) showed that 
the time-resolved spectra have the mean observed a close to — 1 
(i.e. Fe oc E°) and some spectra can be as hard as Fe oc 
E . In spite of these facts, many different versions of t he syn- 
chrotron shock models were proposed r ecently (see e.g. iTavanil 
| l996tlChiang & DermeJl999UPirarll999l) to explain GRB spectra. 
IPanaitescu & Meszaros (2000) hypothesized that inverse Compton 
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scattering of synchrotron self-absorbed radiation can be responsi- 
ble for the observed hard BATSE (i.e. in the range 20-1000 keV) 
spectra under an assumption that electrons emit in a slow cool- 
ing regime, which, however, is h ardly possible in the GRB ejecta 
I Ghisel lini & Celottll999fclGhisellini et all2000h . 

In principle, efficient cooling of electrons can be prevented by 
their reacceleration iLlovd & Petrosiarj200d) . In synchrotron mod- 
els, this requires the fraction of particles taking part in that process 
to be orders of magnitudes smaller than the total number of par- 
ticles (not to exceed the ava ilable energy) and they should always 
be the same iGhisellini et all2000l) . conditions that are difficult to 
imagine. 

Problems with the relativistic synchrotron (and self-Compton) 
models gave rise to optically thick emission models such as quasi- 
thermal Compton i zation IZdziarski & Lambl |l986t iTho mpson 
ll994llLiang||l997HGhisellini & Celottill999l:lsterrJll999l) . where 
the energy is shared among many particles which are now mildly 
relativistic. If the synchrotron radiation is self-absorbed, one can 
achieve rather hard spectra with the peak at 10-50 keV in the co- 
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moving frame of the ejecta. For the bulk Lorentz factor Y ~ 100, 
this peak shifts, however, to an uncomfortably high energy. 

In this paper we show that optically thin synchrotron self- 
Compton mechanism operating at parameters expected in the GRB 
ejecta can naturally produce very hard spectra peaking in the 
BATSE energy band, if the available energy is shared among all 
particles and the particles are continuously accelerated/heated over 
the life-time of a source. In such conditions, the electron/pair distri- 
bution develops a prominent peak at a Lorentz factor 7 ~ 30— 100, 
where synchrotron and Compton losses are balanced by parti- 
cle acceleration and heating due to synchrotron self-absorption. 
A copious pair production reduces available energy per parti- 
cle, moving the spectral peak to lower energies as the burst pro- 
gresses, reproducing t hus the hard-to-soft evolution observed i n 
time-resolved spectra <Ford et aljll995l iRvde & Svenssonl f2002). 
High energ y emission observed in some GRBs (e.g. GRB 941017, 
iGonzalez et alj|2003l). and the prom pt optical emission observed 
in GRB 990123 1 Akerlof et al. 1999) also can be explained in this 
model simultaneously. 



2 MAIN PARAMETERS 

Let us consider the ejecta moving with Lorentz factor Y at the dis- 
tance R from the source. The main parameters determining radia- 
tion physics are the (comoving) energy dissipation rate, magnetic 
field strength B, size of the emission region R', and the number of 
particles described by the Thomson optical depth tt = n c R Vt. It 
is suitable to describe the available energy by the comoving com- 
pactness 

A = -^R'<tt, (1) 

where U is the comoving density of a relevant kind of energy. For- 
mally, A is just the optical depth tt of pairs if we would spend all 
available energy for their mass. In reality, most energy goes to ra- 
diation and the resulting r T < A (e.g., r T ~ 10 - 20 at A = 10 3 
and tt ~ 1 at A = 30). The role of magnetic fields can be de- 
scribed by the magnetic compactness Ab given by Eq. Q with 
U = B 2 /{Six). The ratio Ab/A is model dependent. We assume 
that the magnetic field is below equipartition, i.e. Ab < A. Even in 
magnetically dominating models, one does not expect necessarily 
Ab 2> A, since reconnection of magnetic field providing the en- 
ergy dissipation can reduce its strength within the emission region 
to A B ~ A. 

It is evident that R' should not exceed the size of the causally 
connected region, i.e. R/Y in both transversal and radial direction 
(the latter in the observer's frame is R/Y 2 ). We assume R' = R/Y 
and measure the comoving time t' in units of the light crossing time 
of the region R' /c = R/(cY) which corresponds to the observer 
time R/(cY 2 ). The dissipation compactness (corresponding to the 
energy dissipation rate) is then t — dA/dt' . (For a constant dissi- 
pation rate during t' — 1 we get I — A.) It can be estimated from 
the observed isotropic energy release i? ra d assuming that the en- 
ergy was dissipated homogeneously in a causally connected region 



e = 



m c c 2 Y4nR 2 



7 -E ra d,52r 2 R 1 



(2) 



(Here we adopt notation Q = 10 X Q X in cgs units if not mentioned 
otherwise.) Figure Q demonstrates the levels of the compactness 
and the observed time-scale on a R — Y plane. The observed emis- 
sion episode at the first glance should be a single pulse of time-scale 
R/ (cY 2 ). Actually, most of GRBs have a complex time structure. 
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Figure 1. The map of T — R parameter space. Dotted lines show con- 
stant compactness I levels for energy release -E ra d = 10 52 erg in the 
causally connected region given by Eq. |2}- Dashed lines show the time- 
scale in the observer frame defined as to = R/(cT 2 ). The deceleration 
limit R c i for different £?k is given for the wind with the mass loss rate 
10 — 5 Mq y~ 1 and velocity of 10 3 kms -1 as well as for the interstellar 
medium of nigM = 1 (see Eq. (3))- The circle corresponds to the model 
run simulated in the paper. 



We can then prescribe this episode to a single GRB pulse and asso- 
ciate iJrad with its energy fluence, or to admit that the energy can 
be released by compact flares within the causally connected region. 
In that case, _E ra d should be referred to the fluence of a complex 
emission episode and Eq. then gives an average compactness in 
the region while the local values of I can be substantially higher. 

What is the typical size R where the energy dissipation takes 
place? For radii smaller than < 10 14 cm, the compactness is large 
t > 300 (see Fig.0, pair production is extremely efficient, and the 
mean energy available per particle is rather small. The main emis- 
sion mechanism is then multiple Comp ton scattering. The s pectra 
expected from this regime can be hard tG hisellini & Celottil 19991 
ISternl l 1999), but peaked at too high energies (10—50 keV in the 
comoving frame). At R > 10 cm, the compactness and optical 
depth are smaller, the mean energy available per particle larger and 
the energy is radiated away by optically thin or synchrotron or syn- 
chrotron self-Compton emission. 

The dissipation radius is also limited from above by decelera- 
tion of ejecta in the external environment 



Rd 



1.7 x 10 ie Ekm W3MZ5Y2 2 cm, for wind, 



2.5 x 10" E KB4 n ISM F 2 



2/3 



cm, for ISM, 



(3) 



where the ejecta sweeps up mass ~ M e j/r (see solid lines in 
Fig. 0. Here Ek = YM c jC 2 is the isotropic kinetic energy of 
the ejecta, nisM cm -3 is particle concentration in the interstellar 
medium, Ms is the mass loss rate by the progenitor star in units 
10~ 5 M Q y" 1 , and W3 is the wind velocity in units 1000 kms 
We notice that the deceleration limit in the wind case restricts the 
low compactness (£ < 10) regime to a relatively narrow region of 
parameters. If the kinetic energy is < 10 53 erg, or the environment 
is denser, then the low-compactness regime is hardly possible at 
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all. Thus we consider the typical dissipation radius R — 10 15 cm. 
The observed pulse duration of about 1 s is consistent with the 
fact that GRBs have very l ittle power at time-scales below 1 s 
I Beloborod ov. Stern & Sven sson 2000). 

The optical depth tt in the emission region is bounded from 
above by the opacity of the ejecta 

T oj = 0.3 E K , 54 R^ T^ 1 (4) 

(assuming matter dominated ejecta and that the matter is concen- 
trated within a causally connected shell). For external shocks, tt 
cannot be smaller than that determined by collected material 

_ ( 2 x 10~ 4 Ms W3 1 , for wind, 
Tcxt ~~ \ 2 x 10~ 8 R 17 msM, for ISM. > 

What is the physical mechanism of the energy dissipat ion? I n 
the external shock models iRees & Meszaroslll992t |Pirart [l999). 
the dissipation occurs at R ~ Rd which can be about 10 15 cm 
for the wind environment if Ek ~ 10 53 erg and F ~ 200 (see 
Fig. 0. The impulsive first-orde r Fermi acc eleration would oper- 
ate in a fast cooling regime (see lStenj2003l for a specific version 
of this scenario) which, as discussed in the Introduction, contra- 
dicts the data. An alternative version of shock energy dissipation is 
particle heating by plasma instabilities behind the shock front (see 
e.g. lFrederiksen et all2003l) . We follow this route here. It can also 
operate in internal or 'refreshed' shocks which could be produced 
by collisions of fresh ejecta with previously ejected, partially de- 
celerated material. In the Poynting flux dominated models I Usov 
ll994lLvutikov & Blandford 2003), the magnetic field energy can 
be dissipated at the required distances. 



3 PHYSICAL MODEL 
3.1 Model setup 

As we have seen a number of models can satisfy our requirements. 
Thus we do not specify the exact model for the energy dissipation, 
but consider a toy model where energy is injected to the emission 
region with the constant rate during comoving time R' / c. We adopt 
a slab geometry of the emission region which, in a zero approxima- 
tion, is consistent with both internal and external shock scenarios. 
Indeed, we can expect that the main energy release takes place be- 
hind the shock front in a layer which is thin relative to the size of 
the causally connected region. This geometry can also be a satis- 
factory approximation for the magnetic reconnection scenario. In- 
deed, the magnetic field is probably predominantly transversal and 
the reconnection plane is again perpendicular to the direction of 
propagation. 

We neglect the curvature of the shock front and the bulk ve- 
locity gradient thus reducing the problem to a static slab. In this 
approach we omit a number of effects associated with relativis- 
tic expansion of the emitting shell. These effects are important for 
the description of the time evolution, however they are not critical 
for understanding general spectral properties. The thickness of the 
emitting slab, A < 1 (in units of R'), is unknown since it depends 
on the relaxation process behind the shock front and is probably 
smaller in the case of the reconnection scenario. We take A = 0.1, 
but the results are not very sensitive to its value. 

The energy release is uniform over the slab volume and we 
assume that the energy is injected in a form of acceleration of elec- 
trons and pairs which obtain equal amount of energy per unit time. 
The optical depth can increase due to pair production. We treat the 
magnetic field geometry as chaotic, therefore all pairs are isotropic. 
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The model is fully described by four parameters: (i) the initial 
Thomson optical depth across the slab, To = jieOrA 7?'; (ii) co- 
moving size R'\ (iii) dissipation compactness i\ and (iv) magnetic 
compactness Ab- 

3.2 Radiative processes 

Let us first consider how particles (electron and positrons) of 
Lorentz factor 7 are heated and how do they cool. The energy gain 
rate of a particle is simply given by the heating rate oc i divided by 
the number of particles which is proportional to the total Thomson 
optical depth (including pairs) across the slab tt- Particles cool by 
emitting synchrotron radiation and by scattering this radiation (self- 
Compton mechanism). The energy balance equation can be written 
as: 

S = --^ A b+At)7 2 . (6) 
at tt 6 

Here 77 < 1 accounts for the reduced synchrotron cooling due 
to synchrotron self-absorption, and At is the compactness corre- 
sponding to the energy density of soft photons in the Thomson 
regime (with energy e = E/m c c 2 < I/7). The typical cooling 
time is then t coo \ ~ 1/[(At + ?7Ab)7], which, for the GRB condi- 
tions, is orders of magnitude smaller than the light-crossing time. 
The balance between heating and cooling is achieved at 

7b w v'<7(At + 7 ? Ab)t- 1/2 , (7) 

where At and 77 also depend on 7b. Particles with 7 > 7b 
lose energy faster than they gain it, while at 7 < 7b the situa- 
tion is opposite. As a result, a very narrow electron distribution 
peaked at 7b develops. This allows us to adopt the approxima- 
tion that all particles have the same Lorentz factor 7 = 7b. The 
radiation compactness At can be expressed as a sum of the syn- 
chrotron A s = yryAs and first Compton scattering A c = J/^As 
energy densities (further scattering orders are in the Klein-Nishina 
limit). In the adopted approximation, the Compton parameter is just 
y = £tt7 2 , where the geometrical factor f; ~ 1 for a spherical 
source and £ ~ (2/3) ln(3/2A) ~ 1.8 for a slab with A = 0.1. 
Thus Eq. is reduced to 

y(l + y + y 2 )^Z£/( V A B ). (8) 

When synchrotron self-absorption is negligible, rj = 1, we find the 
solution y w (^/A B ) 1/3 (or y « (^/A B ) 1/2 , if the second 
Compton scattering operates close to the Klein-Nishina limit). At 
small rj, Compton parameter increases. 

The importance of self-absorption depends crucially on 7. The 
optical depth at frequencies below the sync hrotron emission peak 
is (eq. 2.18a in lGhisellini & Svenssorll99ll) 

r s = 15r T /(&7 5 z 5/3 ), a = e/(3 7 2 fo), (9) 

where b = B/B QED and B QED = 4.4 x 10 13 G. Thus the emis- 
sion will be significantly reduced if the self-absorption frequency 
(where r s = 1) is above the emission peak (x > 1). This hap- 
pens at 7 < 7„ = 50(tt,-3/-B3) 1//5 . The same condition for 
tt expressed via Compton parameter is tt > T CI (y) = 0.5 x 

10- 3 (y/O 5/7 B 3 2/7 - 

We now can predict the temporal evolution of the radiation 
spectrum. The optical depth starts growing after about two light- 
crossing times 2A required to produce high energy photons. If 
tt < T CI (yo), the synchrotron is not absorbed and y = j/o, 
7 2 — 2/o/(£tt)> and the synchrotron peak energy decreases with 
optical depth as e s ~ 3j 2 b ~ 2 x 10" 7 (yo/?) 2/7 B\ n t ci /t t ■ 
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The first Compton peak evolves even faster e c i ~ (4/3)7 2 e s ~ 
4 x W- 4 (yo/0 4/7 B'i /7 (T cv /T T ) 2 . When r T grows above r cr , r/ 
decreases because of self-absorption, y correspondingly increases 
(Eq. [8j, the synchrotron peak becomes more stable, and the re- 
sulting electron energy and the first Compton peak start to evolve 
slower. The second Compton peak is produced in the Klein-Nishina 
limit at small tt, and evolves slowly e C 2 ~ 7 oc 1/^/tt, while at 
larger t t the evolution speeds up. We now check there predictions 
by numerical simulations. 



3.3 Simulations 

The simulations were performed u sing a Large Particle Monte- 
Carlo code (LPMC) developed by ISternl <1985l) and IStern et alJ 
1 1995). It treats Compton scattering, photon-photon pair produc- 
tion and pair annihilation, synchrotron radiation and synchrotron 
self-absorption. The exact cross-section are used for the first three 
processes, while the cross-sections in the relativistic approxima- 
tion are used for the synchrotron process iGhisellini & Svenssonl 
Il99ll) . The electron/pair and photon distributions are computed 
self-consistently. The code is essentially nonlinear: the simulated 
particles constitute at the same time a target medium for other par- 
ticles. The geometry of the emission region is a pill-box of radius 
R' and thickness A = 0.1. Output photons are recorded when they 
cross the surfaces z = ±(A/2 + 0.1)/?'. 

As an example, we take typical parameters described in § [2] 
The comoving radius R' — 10 13 cm, the initial Thomson optical 
depth to = 6 X 10~ 4 (close to the critical r cr , but higher than r cxt ), 
compactness £ — 3 (corresponding to T « 130 for £^,52 = L 
see Eq.[2}, and the magnetic compactness Ab = 0.3 (correspond- 
ing to the comoving magnetic field B' ~ 1000 G, three times be- 
low the equipartition). The evolution of broad-band spectra and the 
electron distribution are shown in Fig. [2] 

At the start of simulations, t t = to < r cr and the electron 
heating and cooling are balanced at high 7 (see dashed curves). 
Partially self-absorbed synchrotron radiation (solid curves, lower 
energy bump) peaks at e s ~ 3 x 10 ~ 7 in the comoving frame an d 
has a hard low-energy tailF B oc £ 2 lGhisellini& Svensson 1991). 
The optical depth grows nearly linearly with time due to electron- 
positron pair production (see Fig.|3t) and the mean particle energy 
decreases as (7) oc \j\pT (Fig-[3j)). The Compton parameter y = 
£ t t(7 2 ) computed from the pair distribution (see Fig.[3J:) is about 
8 in the beginning, reaches minimum at t' = 0.3 and grows to 
more than 10 by the end of energy injection. The 'observed' y- 
parameter, i.e. the ratio of luminosities in the first Compton bump 
to the synchrotron component grows monotonically, however, from 
3 at t' = 0.1 to 16 at t' = 1. This discrepancy is caused by the non- 
stationarity of the problem - it takes too long time to build up the 
spectrum and to reach a steady-state. 

The first Compton peak oc (7 2 ) moves to softer energies and 
crosses the 'BATSE window' (note that due to self-absorption, syn- 
chrotron peak energy is very stable). This spectral evolu tion is con- 
sistent with the observed in t ime-resolved pulses (e.g. iFord et all 
ll995tlRvde & SvenssonlioO^ . The photon flux in the BATSE band 
(Fig. |5Ji) shows a 'fast rise - exponential decay' behaviour often 
seen in GRBs. It decays much before the energy supply terminates, 
but has a long flat part. While the second Compton peak at ~ 10- 
100 MeV, on the contrary, rises later and decays on longer time- 
scale. 

We fitted the model photon spectrum Ne = Fe/E in the 
BATSE window by a phenomenological 'GRB function' consisting 
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Figure 2. Instantaneous (comoving frame) photon spectra within the slab 
(i.e. the source functions) and corresponding electron distributions. Solid 
curves show the photon spectra EFe (in arbitrary units) at times of 
0.1,0.3, 0.6, and 1 (in units R'/c). The spectra consist of a low-energy 
(partially self-absorbed) synchrotron bump and two Compton scattering or- 
ders. Further scatterings are suppressed due to the Klein-Nishina effect. 
The electron energy distribution function 7(7 — l)d7V/d7 for the same 
time intervals is shown by dashed curves (the peak evolves towards lower 
energies). Parameters of simulations: R' = 10 13 cm, to = 6 X 10~~ 4 , 
1 = 3, Ab = 0.3. Dotted line shows the hardest possible Fe oc E 1 
power-law reachable at the low-energy slope of the first Compton bump. 
The BATSE (20-1000 keV) and EGRET (3-100 MeV) bands redshifted to a 
comoving frame by (1 + z)/2V (where 2 ~ 1 is the cosmological redshift 
and r = 130) are shown by horizontal bars. 

of a power-law with an exponential cutoff, Ne oc E a e xp(— E(2+ 
a)/E n ), merging to a high-energy power-law oc E 13 I Ba nd et all 
1993). The evolution of the fitted (observed) peak energy E p of the 
EFe spectrum and the spectral slopes is shown in Fig.|3^,f. One 
sees softening of the spectrum as the burst progresses. The fitted a 
is close to —1, the most probable value in the distribution of time- 
resolved spectra jPreece et all200d) . The results of spectral fitting 
depend somewhat on the used energy interval: a becomes larger 
(spectrum hardens) in a wider interval and /3 is softer. Often the 
data at lower energies (with better statistics) dominate the fitting 
procedure, then the fitted a can be much harder. A correlation be- 
tween a and E p is also expected. The hardest possible spectrum, 
q = 0, correspo nds to Compton scattering by isotropic mo noen- 
ergeti c electrons iBlumenthal & Gouloll97ol : lRvbicki & Lightmanl 



4 SUMMARY 

The radiative processes that are responsible for particle cooling de- 
pend on the optical depth in the emission region. If tt is very small, 
the mean 7 is large and the BATSE photons would be produced by 
synchrotron radiation. At tt ~ 1, quasi-thermal Comptonization 
operates. At intermediate tt ~ 3 x 10~ 4 -10~ 2 , the electrons/pairs 
have 7 ~ 30 — 100 and the first Compton peak is produced in 
the BATSE band. Such optical depth seems natural for the external 
shock in the typical Wolf-Rayet progenitor wind as well as for the 
emission within the ejecta (e.g. due to magnetic reconnection or 
collisions of the fresh ejecta with already decelerated material). Ef- 
ficient pair production at intermediate compactnesses, I = 0.1-10, 
also can be responsible for the required tt. 

Synchrotron self-Compton emission from continuously 
heated, nearly monoenergetic electrons can explain many observed 
features of GRBs. The synchrotron spectrum peaking at e ~ 
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Figure 3. Evolution (a) of the Thomson optical depth due to pair produc- 
tion, (b) of the mean electron Lorentz factor (7), and (c) of the Compton y- 
parameter. (d) Photon flux in the BATSE band (solid curve) and the EGRET 
band (dotted curve) in arbitrary units. Evolution of spectral parameters (e) 
E p , (f) a (circles), and /3 (squares). 



3 x 10 7 in the comoving frame will be blueshifted to the ex- 
treme UV region. The self-absorbed low-energy tail is hard, Fe oc 
E 2 , and can explain the p rompt optical radiation detected from 
GRB 990123 lAkerlof et ailll999l) . The energy of the first Comp- 
ton peak is expected to decrease as the burst progresses since pair 
production reduces the mean energy available per particle. In the 
case of a larger To and/or a higher compactness and/or a smaller V, 
the first Compton component peaks in X-rays and possibly can be 
identified with the observed X-ray flashes (e.g. Heise 2003). 

The GRB spectral hardness distribution IPreece et alj f2000) 
can also be reprodu ced (maybe except its hardest events, see 
lGhirlandaet"aill2 003). Since the incident synchrotron spectrum is 
hard, the low energy slope of the scattered radiation, Fe oc E , 
is determined by kinematics of single Compton scattering by mo- 
noenergetic electrons. 

The second inverse Compton peak observed at 10 MeV - 10 
GeV is delayed relative to the soft gamma-ray emission and lasts 
longer. In spite of large y-parameter, it does not necessarily dom- 
inate the total energy output because of the Klein-Nishina effect. 
The rather hard (a ~ — 1) spectrum at ~ 10 MeV can match ob- 
servations of GRB 941017 <Gonzalez et all2003h and its observed 
slow evolution. If the mean particle energy decreases rapidly, this 



component can possibly produce even the second distinct pulse in 
the BATSE range. 

Summarizing, the proposed model can explain a large fraction 
of GRB spectra and their time evolution. It also reproduces the high 
energy ~ 10 — 100 MeV emission detected by EGRET in some 
bursts and the prompt optical emission. Two latter phenomena are 
natural within this model and do not require additional assumptions 
or separate emission regions. 
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